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ABSTRACT 

We present for the first time phase-resolved UV spectroscopy of an Of?p star, 
namely, HD 191612. The observations were acquired with the Space Telescope Imag- 
ing Spectrograph (STIS) on-board the Hubble Space Telescope (HST). We report the 
variability observed in the main photosphcric and wind features and compare the re- 
sults with previous findings for the Of?p star HD 108. We show that UV line strengths, 
Ha, and longitudinal magnetic field, vary coherently according to the rotational pe- 
riod (P ro t = 537. 6d), providing additional support for the magnetic oblique rotator 
scenario. The stellar and wind parameters of HD 191612 are obtained based on NLTE 
expanding atmosphere models. The peculiar wind line profile variations revealed by 
the new STIS data - not reproduced by ID atmosphere models - are addressed through 
non-spherical MHD simulations coupled with radiative transfer. The basic aspects of 
the UV variability observed are explained and the structure of the dynamical magne- 
to sphere of HD 191612 is discussed. 

Key words: stars: winds - stars: atmospheres - stars: massive - stars: magnetic 
fields. 



1 INTRODUCTION 

The peculiar Of?p stars have been studied in detail in the 
optical part of the spectrum, with observations spanning 
several decades for HD 108 and HD 191612 (e.g., Walborn 
1973; Naze et al. 2010). Large-amplitude variability of op- 
tical features (e.g., of Ha and He I lines) is found to be 
common for this spectral class (Naze, Vreux & Rauw 2001; 
Howarth et al. 2007), and is characterized by 'high' and 
'low' (or maximum/minimum) states in the degree of emis- 
sion. This behaviour is now understood in the context of a 
magnetic oblique rotator model, wherein the observer's view 
of a magnetically constrained stellar-wind outflow is modu- 
lated by stellar rotation (see Donati et al. 2006; Wade et al. 
2011; Sundqvist et al. 2012). 
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Despite advances in understanding the nature of Of ?p 
stars, the behavior of the ultraviolet (UV) spectra at differ- 
ent optical states remains poorly known. Such information is 
crucially important, since this spectral region contains the 
most sensitive diagnostics of hot-star winds. We therefore 
lack, so far, a complete view of the stellar-wind confinement 
phenomenon (ud-Doula & Owocki 2002) in these stars. A 
first effort to address this problem was recently reported by 
our group for the very slow rotator HD 108 (Marcolino et 
al. 2012; hereinafter Paper I). New observations in the ultra- 
violet were acquired close to the optical minimum state; the 
UV spectrum did show changes compared to International 
Ultraviolet Explorer (IUE) observations acquired during a 
high state, but the variability was not as dramatic as ob- 
served in the optical. Marcolino et al. (2012) concluded that 
the stellar wind of HD 108 is not severely affected by the 
magnetic field on a global scale (i.e., at large radii; see Pa- 
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per I for details). However, because of the long rotational 
period of HD 108 (probably about 55 years; Naze, Vreux & 
Rauw 2001; Naze et al. 2010), we still lacked information 
about UV line profiles at various rotational phases. Thus, 
until now, we have not had a detailed knowledge of UV 
variability matching that in the optical, for any Of?p star. 

In this paper we present high-resolution UV spec- 
troscopy of HD 1916133 (Piot = 537.6 d; Howarth et al. 
2007), which represents the first phase-resolved UV spectra 
of an Of?p star covering the entire rotational cycle. The ob- 
servations are described in Section [21 where we also summa- 
rize the variability observed in the main photospheric and 
wind features and compare the results with our previous 
findings for HD 108. We demonstrate that the strengths of 
the main UV wind lines, Ha, and longitudinal magnetic field 
(Bz), vary coherently according to the rotational period, 
thereby providing support for an origin linked to the mag- 
netic oblique rotator model. The stellar and wind parame- 
ters of HD 191612 are then revised and discussed based on 
UV modelling using non-LTE expanding-atmosphere mod- 
els (Section [3]). The peculiar variability revealed by the 
new STIS data is not reproduced by ID models, but can 
be explained qualitatively by using magnetohydrodynamical 
(MHD) simulations of a magnetically confined wind (Section 
[4| . Section [5] summarizes the main results of our work and 
discuss their consequences for the Of?p class as a whole. 



2 OBSERVATIONS 

We used the Space Telescope Imaging Spectrograph (STIS) 
on board HST to obtain ultraviolet spectra of HD 191612 
(PI. Bouret; program ID 12179). Table [1] provides a log of 
the observations, which were obtained at four phases sepa- 
rated by a quarter of the star's 537.6d rotational cycle. At 
each phase, the observations consisted of two pairs of ex- 
posures through the small (0.2" x 0.2") aperture, first with 
the E140M echelle grating and then with the E230M echelle 
grating. 

Observations with the E140M grating were centered at 
1425 A and recorded by the FUV MAMA detector, while ob- 
servations with the E230M grating were centered at 2707 A 
and recorded by the NUV MAMA detector. The total wave- 
length coverage is ~1150 - 1750 A (E140M) and ~2250 - 
3100 A (E230M), with resolving powers R = 45,800 and 
30,000, respectively. The modified Julian Dates (MJD) of 
the mid-exposures are listed in Table [T] along with the to- 
tal integration time (divided equally between the two expo- 
sures). Since the two exposures obtained with each config- 
uration were combined, the mid-point of the sequence was 
used to calculate the phase (j> according to the ephemeris 
published by Howarth et al. (2007). Due to the very long 
stellar variability period, the NUV/FUV exposure pairs cor- 
respond to essentially identical phases (to within a few 10 -5 
cycles). 



1 HD 191612 is a known double-lined spectroscopic binary, with 
P orb = 1542d (Howarth et al. 2007). 



Table 1. Log of HST/STIS observations of HD 191612. Rota- 
tional phases are from the ephemeris given by Howarth et al. 
(2007): JD = 2453415.2 ± 537.6B. The phase uncertainty is about 
0.002. The total exposure time is given in seconds. 



Grating 


Date (UT) 


MJD (mid) 


exp. time 


Phase tj> 


E140M 


2010-08-23 


55431.089 


1500 


0.751 


E230M 


2010-08-23 


55431.104 


300 


0.751 


E140M 


2011-01-04 


55565.294 


1500 


0.000 


E230M 


2011-01-04 


55565.310 


300 


0.000 


E140M 


2011-05-18 


55699.291 


1500 


0.250 


E230M 


2011-05-18 


55699.303 


300 


0.250 


E140M 


2011-09-30 


55834.512 


1500 


0.501 


E230M 


2011-09-30 


55834.527 


300 


0.501 



2.1 Observed variability characteristics 

In Figure[T]we present the spectra of HD 191612 obtained at 
different rotational phases. In this section we focus mainly 
on FUV data, where we find the main stellar wind diagnos- 
tics (N v A1240, Si iv A1400, C iv A1550, and N iv A1718) 
and the iron "forest" . Data obtained in the NUV region re- 
veal fewer lines which are mostly in absorption. The most 
prominent is C in A2296, which we present in Figure [21 

There is significant line-profile variability at a level 
which exceeds that normally found in O-star spectrfl al- 
though the amplitudes are nonetheless modest compared 
to the substantial changes reported in the optical spec- 
trum (see, for example, Fig. 5 of Walborn et al. 2003). The 
main changes are observed in the absorption troughs of the 
N v A1240 and C iv A1550 P Cygni profiles, together with 
the high-velocity part of N iv A1718, which are all stronger 
at <j> = 0.0 (the 'high state') than at other phases, while 
the (essentially photospheric) iron-forest absorption weak- 
ens at this phase (Figure [1]). In contrast, both Si iv A1400 
and C ill A2296 are weaker at <j> = 0.0 (see Figure [2}. The 
C ill A1176 multiplet (not shown) is very weak compared 
to these two features but follows the same trend. This be- 
haviour is related to subtle magnetic field effects on the wind 
and is discussed more in detail in Section [4] 

The Of?p star HD 108 also exhibits greater wind ab- 
sorption in C iv and N v and weaker absorption in the 
iron forest during the high state. However, there are im- 
portant differences between the two stars. For example, in 
HD 191612 (unlike HD 108) the C iv and N v wind pro- 
files never become saturated, regardless of the phase. More- 
over, in HD 191612 (again, unlike HD 108) the Si iv A1400 
behaves in the opposite way to C iv and N v, the absorp- 
tion being strongest at phase 0.5. Such differences may be 
related to the different stellar and magnetic properties of 
these stars. Different mass-loss rates, magnetic obliquity, 

2 Many O stars are known to exhibit resonance-line variability 
in the form of discrete absorption components (DACs); see, e.g., 
Kaper et al. (1996). DACs are localized optical depth enhance- 
ments that propagate through P-Cygni absorption troughs from 
red to blue on a timescale of typically ~days, becoming narrower 
as they approach the terminal velocity. The variations seen here 
are fundamentally different in character, insofar as the change in 
absorption depth occurs over the entire width of the line profile. 
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Figure 1. Phase-resolved data of HD 191612: <f> = 0.0, 0.25, 0.50, and 0.75. The continua are essentially the same, the total energy 
emitted in the UV is essentially preserved. The N V and C IV profile changes from rf> = 0.0 and 0.5 follow the trend observed in HD 108 
(high to low state; Marcolino et al. 2012). Flux units are erg cm -2 s -1 A -1 . Wavelengths are given in angstroms. 
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Figure 2. The C III A2296 feature at different rotational phases 
- NUV/MAMA STIS spectra. Flux units are erg cm -2 s _1 A -1 . 



Table 2. The high-state (<f> ~ 0.0) behaviour of the main UV 
features in HD 191612 and HD 108 compared to low-state spectra. 
The +abs (— abs) notation indicates greater (less) absorption in 
the photospheric profiles or in the absorption part of a P Cygni 
profile in the high state. Note the different behaviours for Si IV 
and C III. 



Feature 


HD 191612 


HD 108 


N V A1240 


+abs 


+abs 


C iv A1550 


+abs 


+abs 


N iv A1718 


4-abs 


+abs 


Fe iv forest 


—abs 


—abs 


C in A1176 


-abs 


? 


Si iv A1400 


-abs 


-|-abs 


C m A2296 


-abs 


? 



and hence viewing angle may all affect the observed line 
profile variations. 

Table [2] summarizes this discussion of the relative vari- 
ability of the main UV features of HD 108 and HD 191612. 
Note that the analysis of the C ill lines in HD 108 is hin- 
dered by the very low signal-to-noise of the IUE observations 
available. 



2.2 Equivalent width measurements 

To quantify the variability, we measured equivalent widths 
for the main UV lines, using fixed wavelength intervals for 
each feature (e.g., 1534.0-1551.3A for C iv and 1390.2- 
1395. OA for Si iv). In Fig. |3]we present the results for C iv 
and Si iv, along with Ha equivalent widths and the longitu- 
dinal magnetic-field values, Bz (from Howarth et al. 2007; 
Wade et al. 2011). Maximum C iv A1550 and minimum 
Si iv A1400 absorptions occur at cj> = 0, which coincides 
with maximum Ha emission. The UV iron-forest strength 
(not shown), as well as other UV profiles (e.g. N v A1240 
and C ill A2296), also follows the rotational variability, along 
with photometry from Hipparcos (Walborn et al. 2004). In 
summary, the UV profiles of HD 191612 vary in concert 
with other measurements: the UV variability follows the es- 
tablished 537. 6-d period and results from rotational modu- 
lation. 

The oblique rotator model was proposed for HD 191612 
by Donati et al. (2006), with details subsequently refined 
by Howarth et al. (2007) and Wade et al. (2011). The cur- 
rently accepted model has a (rotational) axial inclination of 
i ~ 30°, a centred dipole magnetic field with an obliquity 
of j3 ~ 67°, and a polar strength of ~2.5kG. With this ge- 
ometry, the line of sight is almost aligned with the dipole 
axis at rotational phase (f> ~ 0, and almost perpendicular to 
it at ~ 0.5. In such a scenario, we expect enhanced den- 
sity and slower velocities in the region of the magnetic equa- 
tor, which is fed by magnetically deflected outflows from the 
polar and temperate latitudes (ud-Doula & Owocki 2002). 
Thus, a priori, the net effect expected for a wind profile at 
4> — 0.5 is less absorption at high velocities, as exhibited, 
e.g., by C IV A1550 and N v A1240. At a first glance, the 
Si iv A1400 profile (and C iv A2296) seems to contradict 
this idea. We come back to this question later in Section 3J 
where we show that effects from a dynamical magneto sphere 
may solve this puzzle. 



3 ATMOSPHERE MODELS 

For the oblique-rotator configuration described above, 
spherical symmetry is clearly a very poor approximation 
for lines principally formed inside the Alfven radius; it is 
precisely this lack of spherical symmetry (coupled with for- 
tuitous viewing geometry) that generates the strong rota- 
tional modulation observed in Ha and some optical helium 
lines. This conclusion is strengthened by the success that 
Sundqvist et al. (2012) had in reproducing the observed 
variability of Ha by using the structure of a magnetically 
confined wind predicted by MHD models. 

Sundqvist et al. (2012) coupled their MHD models to 
relatively simple radiative-transfer calculations appropriate 
to Ha line formation. However, this work has not yet pro- 
gressed to the point where fully self-consistent modelling of 
all UV lines is possible (see however Section [3]). In order to 
obtain a quantitative characterization of our STIS data, we 
therefore first turn to the other extreme of spectral mod- 
elling, embodied in the CMFGEN code (Hillier & Miller 
1998). This code solves the equation of radiative transfer 
with considerable attention to the microphysics (line blan- 
keting, NLTE, radiative and statistical equilibria), but with 
simplified - normally spherical - symmetry. 

While spherical symmetry is obviously a poor assump- 
tion for lines formed in the vicinity of the Alfven radius, we 
might expect that features formed in or close to the pho- 
tosphere, such as the iron forest, can reasonably be mod- 
elled under this approximation. Furthermore, because the 
magnetic dipole energy density falls much faster with radial 
distance than does the wind kinetic energy density (~ r~ 6 
vs. r~ 2 ), the stellar-wind outflow must become essentially 
radial at large distances. With some caution, we may there- 
fore investigate lines formed at large radii - such as the high- 
velocity parts of UV P Cygni profiles - under the simplifying 
assumption of spherical symmetry. 

In order to further avoid non-spherical effects, we fo- 
cused on the analysis of STIS data at phase = 0, when 
the magnetic pole is oriented towards the observer (see Fig. 
[3]). At this configuration, the line-of-sight absorption column 
will consist mostly of wind plasma streaming from the open 
field lines above the magnetic pole, which presumably then 
leads to a reduced influence of the dense, confined wind ma- 
terial concentrated near the magnetic equator. Thus, for UV 
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Figure 3. Phase-resolved data of HD 191612. Upper panels: Longitudinal magnetic field component (left) and Ho equivalent width 
(right; data from Howarth et al. 2007; Wade et al. 2011). Lower panels: equivalent width of C IV A1550 and Si IV A1400 obtained from 
the new STIS observations. 



wind lines, a spherically symmetric model should be a better 
approximation in the high state. 

In Table [3] we summarize the stellar and wind proper- 
ties inferred. Our final model is compared in Figure [4] with 
normalized STIS spectra. The effects of clumping were ac- 
commodated with a filling-factor formalism (f c i — 0.1), be- 
ginning at an expansion velocity of 30 km s _1 (v c i); we also 
allowed for the the effects of X-rays produced by wind in- 
stabilities, adopting logLx/^Boi = —6.1 ±0.1 (Naze et al. 
2010). 

Below we discuss other parameters and the issues found 
during the analysis. 



3.1 Chemical composition, rotation and 
macroturbulence 

Following Howarth et al. (2007), we adopted a macrotur- 
bulent velocity v mac = 45 km s^ 1 and projected equatorial 
rotation velocity v e smi of zero. We note, however, that the 
UV observations are equally well matched by any model 
having 



\J »L + (« e sin i/2) 2 ~ 45km s 1 . 

We used the solar chemical composition (Grevesse et 
al. 2010), except for nitrogen. The nitrogen content in 
HD 191612 was recently determined by Martins et al. (2011) 
based on the analysis of photospheric lines in high-resolution 
optical data. Here, we use their results, which fix N/H at 
2.7 x 10~ 4 (the solar value is 6.8 x 10 -5 ; by number). 

3.2 Luminosity and distance 

HD 191612 is probably a member of the Cyg OB3 associa- 
tion, at a distance of 2.3 kpc (Humphreys 1978). Adopting 
this distance, we estimate the luminosity from a fit to the 
FUV and NUV spectra plus UBVJHK magnitudes available 
in the Galactic O-Star Catalog (Sota et al. 2008). The re- 
sults are included in Table [3] and the fit displayed in Fig. [S] 
The differential extinction E(B — V) was estimated from the 
observed B — V — 0.26 and the typical intrinsic value for 
06-08 stars, (B - V)o = -0.27 (Martins & Plez 2006). 

We note that HD 191612 is known to be photometri- 
cally variable, but the amplitude is too small to affect these 
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Figure 4. Our adopted CMFGEN model (red) compared to STIS data of HD 191612 at phase <fi = 0.0 (black). Wavelengths are given 
in angstroms and the flux is normalized. H I interstellar absorption is accounted for shortward the N V profile; log N(H i) = 21.48. 




Figure 5. STIS FUV and NUV data plus UBVJHK photometry converted to flux points (black) and our adopted CMFGEN model 
(red). The adopted distance is 2.3 kpc. Flux units are erg cm -2 s _1 A -1 . 
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Table 3. Revised stellar and wind properties of HD 191612 (CM- 
FGEN models). 



T .// (K) 
log g (cgs) 
log L*/L 
R*/R© 
M*/M 

v e q sini (km s -1 ) 
(km s" 1 ) 



cphot 



Prot (days) 



(km s 1 ) 



36000 ± 2000 
3.75 ± 0.20 
5.45 ± 0.10 
i q 7 +2.3 

10. I _2 

39 tf 7 


45 
15 

537.6 



log M 


-7.9 ± 0.3 


Voo (km s _1 ) 


2400 ± 500 


f c ( 


0.1 


v ci (km s _1 ) 


30 


log Lx/Lbol 


-6.1 


distance (pc) 


2300 



more than ~0.05. We also neglect the contribution of the 
~early-B secondary. Such a star is much cooler (~20kK) 
and fainter than the primary and is expected to contribute 
only about 10 per cent in the visible, and considerably less 
at UV wavelengths (see Howarth et al. 2007). 



3.3 Issues 

From Fig. [4] we see that the agreement between our atmo- 
sphere model and the STIS data of HD 191612 at phase 
(j> = is reasonable. We predict P Cygni profiles for 
N v A1240 and C IV A1550, a good fit for the iron forest, 
and absorption/incipient P Cygni profiles to Si IV A 1400 and 
N IV A 1718, as observed. However, it is evident that some 
features could not be reproduced in detail. For example, the 
synthetic N IV A1718 profile lacks blueshifted absorption. By 
increasing the N/H ratio to ~ 8 x 10 -4 we improve the fit 
to the observed profile, but this value is much higher than 
that reported by Martins et al. (2012), even considering the 
uncertainties in the determination. We have also tried sev- 
eral ways to improve the fit to the Si IV A1400 feature (see 
Fig. g}, although a comparison is handicapped by interstel- 
lar absorption. 

Some discrepancies were already expected in our analy- 
sis, since we used spherically symmetric atmosphere models. 
As discussed previously, lines formed closer to the Alfven 
radius may be significantly affected by non-sphericity. Ide- 
ally, MHD wind-confinement simulations coupled with de- 
tailed 3D radiative-transfer calculations should be used to 
derive the stellar and wind properties of HD 191612. How- 
ever, given the complexity of such a task, the calculation of 
a full UV and optical theoretical spectrum - as is done by 
ID models - is not currently feasible. Instead, in the next 
section we attempt to provide a physical interpretation for 
generic UV profiles based on MHD simulations. 



4 UV WIND LINES CALCULATED FROM 
MHD MODELS 

This section reports initial results from ongoing attempts 
to use radiation MHD wind simulations (developed first by 
ud-Doula & Owocki 2002) to analyze the strengths and vari- 
ability of UV resonance lines in magnetic massive stars. To 
this end, we apply the same MHD simulation used previously 
to model the rotational phase variation of HD 191612's Ha 
line (Sundqvist et al. 2012), however modifying the radia- 
tive transfer code developed there to synthesize UV scatter- 
ing lines (following the '3-D SEP method, e.g. Cranmer & 
Owocki 1996). 

Fig. [5] shows generic synthetic profiles for a weak (up- 
per panel) and moderately strong (lower panel) line, calcu- 
lated i) from the MHD model for an observer viewing from 
above the magnetic pole ("high state") and equator ("low 
state"), and ii) from the spherically symmetric wind model 
used as initial condition to the MHD simulation. Note how 
the phase trends of the modeled absorption-strengths indeed 
display opposite behaviour for the weak and strong lines, as 
observed for the (weaker) Si IV and (stronger) C IV lines in 
HD 191612 (see Fig. [[J. 

This anti-correlated variability of strong and weak lines 
can be qualitatively understood as follows: the line-of-sight 
absorption column in the low state samples the low-velocity, 
high density wind material trapped by closed field lines near 
the magnetic equator in HD 191612's dynamical magneto- 
sphere (Sundqvist et al. 2012). By contrast, the absorption 
column in the high state covers a wider range of velocities 
and densities, since it stems from the open field lines stream- 
ing from the magnetic pole, and so more closely resembles 
a 'normal' wind outflow. The result is that a strong line, 
which has sufficient optical depth to produce measurable 
absorption at high velocities, is stronger in the high state 
(lower panel). On the other hand, a weak line, which is not 
able to produce significant absorption at high velocities, is 
stronger in the low state (upper panel) , because during that 
phase the line can absorb from the high density, low velocity 
material trapped in the wind's dynamical magnetosphere. 

This figure also shows that profiles computed for in- 
trinsically strong lines (e.g., those corresponding to the line 
strength parameter «o = 1 in Fig. [6)l are systematically 
weaker than their counterparts computed from spherically 
symmetric models, particularly at high velocities. 

This effect can be understood by considering the faster- 
than-radial expansion of a radial polar flow-tube in a mag- 
netic wind (Kopp & Holzer 1976; Owocki & ud-Doula 2004), 
which leads to a different density profile as function of veloc- 
ity than in a spherical wind. Following Owocki & ud-Doula 
(2004), let us parameterize the non-radial expansion factor 
of an areal element, A(r'): 



h(r') 



A(r') _ r'(l + v l /s ) 



r' + rj 



3/8 



where r' = r/R+ is the scaled radius and 77* is the mag- 
netic wind confinement parameter applicable to a dipolar 
field (ud-Doula & Owocki 2002). Note how h starts off at 
unity, but increases with radius reaching the asymptotic 
value h(r' —¥ 00) — 1 + rj^ S . By the equation of continuity 
p ~ l/(Av); consequently, for equal basal mass fluxes, the 
density p of a non-magnetic wind will be higher at a speci- 
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Figure 6. Synthetic UV scattering lines computed from MHD 
wind simulations for an observer viewing the star from above 
the magnetic pole (blue, dashed lines) and equator (red, dashed- 
dotted lines), as compared to profiles computed from a spherically 
symmetric, non-magnetic wind model (black lines). The lower 
(upper; shifted by +1.5) panel shows profiles computed using a 
constant line-strength parameter kq = 1 [kq = 0.1), where reo is 
proportional to the wind mass-loss rate and the ion fraction (Puis, 
Vink, & Najarro 2008). The abscissa shows the wind velocity 
normalized to the terminal speed of the non-magnetic model. 



fied velocity in the outer wind than the density pe of a mag- 
netic one, by a factor p/pB ~ (1 + V*^ 8 )- Since r] t ~ 50 for 
HD 191612, p/pB ~ 5.3. This effect may partly explain the 
low mass-loss rate inferred from CMFGEN models (Sect. 
dilution effect neglected) compared to the rate needed in the 
MHD simulations to reproduce HD 191612's observed Ha 
emission strength and variability (~ 10 -6 Mq yr _1 ). Note 
also from Fig.[6]that the MHD simulation further predicts a 
significant increase in wind terminal speed, here by ~ 50%. 
However, this does not seem to be observed for HD 191612, 
whose Voo ~ 2400 km/s is not much higher than typical val- 
ues inferred for non-magnetic O-stars. 

A detailed study of the formation of UV wind lines in 
magnetic massive stars, focusing on issues like those dis- 
cussed here, will be reported separately (Sundqvist & ud- 
Doula, in prep.). 



5 SUMMARY AND CONCLUSIONS 

We have presented for the first time phase-resolved UV data 
for an Of?p star, HD 191612. Our main findings and con- 
clusions are outlined below. 

• The spectra observed at different rotational phases are 
broadly similar, with the same set of wind and photospheric 
lines consistently present. However, line profile variations 
are very significant compared to the ones found in normal 
O stars. 

• The variability observed is not as drastic as that re- 
ported for Ha and the He I lines in the optical part of 
the spectrum. We interpret this result as indicative of the 
stronger influence of the magnetic field on the lines formed 



closer to the photosphere (i.e., inside the Alfven radius), 
where we find confined material. 

• We observe more wind absorption in N v A1240 and 
C IV A1550 and a less intense iron-forest spectrum at <f> = 
compared to cj) = 0.5. The former phase corresponds ap- 
proximately to the configuration where the magnetic pole 
is towards the observer while the latter corresponds to the 
line of sight lying in the plane of the magnetic equator. Such 
variability is very similar to that reported for HD 108 in Pa- 
per I, and could be considered as probably representative of 
Of?p stars. There are, however, notable differences between 
HD 108 and HD 191612, such as the absence of saturated 
P Cygni profiles in HD 191612 and the anti-correlated vari- 
ability of Si IV A1400 and other wind lines. 

• CMFGEN models are capable of providing a reasonable 
fit to the UV spectrum of HD 191612 at phase <f> = 0. The 
wind and photospheric lines are correctly predicted in broad 
terms, but not all profiles are reproduced in detail. Further- 
more, spherically symmetric models cannot reproduce the 
variability revealed by the STIS data, since these variations 
arise from the dynamical magnetosphere of HD 191612. 

• We use MHD simulations in combination with simpli- 
fied radiative transfer calculations to analyse the new STIS 
phase resolved data of HD 191612. We compare generic UV 
profiles representing C IV A1550 and Si IV A1400. The ob- 
served behaviour at phases cj> — 0.0 and <f> = 0.5 are qualita- 
tively reproduced for both lines; an explanation for the anti- 
correlated variability (Figs. [T] and is provided in terms of 
line strengths. 

• The mass-loss rate derived for HD 191612 from non- 
LTE expanding atmosphere models for the UV is consider- 
ably lower than the one inferred from the MHD simulations 
(based on Ha; ~ 10 -6 Mq yr _1 ). This discrepancy can be 
partially understood by a density dilution effect (Sect. [4| 
but deserves a deeper investigation. 



The present results take us a further step towards a de- 
tailed picture of the optical and UV variability in Of?p stars. 
In general, the optical spectra are strongly affected by the 
magnetic field, with line formation at or below the Alfven 
radius. The UV P Cygni profiles present moderate changes, 
suggesting spherical symmetry may be a reasonable first ap- 
proximation for the wind at large radii. MHD simulations 
such as the ones computed here are however mandatory to 
describe the (phase-resolved) UV profiles in detail. 

Two additional Of?p stars have been recently detected 
to be magnetic, namely, NGC 1624-2 and CPD -28 2561 
(Walborn et al. 2010; Wade et al. 2012; Hubrig et al. 2012). 
The basic results regarding the ultraviolet spectrum re- 
ported here and in our previous paper are expected to be 
also observable in these stars. However, for NGC 1624-2, fu- 
ture UV observations will likely reveal a considerably larger 
variability in the wind profiles than the ones seen so far, 
since its dipole field is estimated to be about 20kG (Wade 
et al. 2012). Such a field is about 8 times more intense than 
in any other known magnetic O-type star. 

The Of?p class now stands out among the O-type stars 
as magnetic variables. They represent therefore an ideal as- 
trophysical laboratory for the study of stellar wind confine- 
ment and magnetic field effects in the evolution of massive 
stars. The progress achieved so far shows that multiwave- 
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length analyses are crucial for a complete description of their 
properties. 
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